Nonequilibrium critical relaxation in the presence of extended defects 
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We study nonequilibrium critical relaxation properties of systems with quenched extended de- 
fects, correlated in dimensions and randomly distributed in the remaining d — Ed dimensions. 
Using a field-theoretic renormalization-group approach, we find the scaling behavior of the nonequi- 
librium response and correlation functions and calculate the initial slip exponents 6 and 9' , which 
describe the growth of correlations during the initial stage of the critical relaxation, in the two-loop 
approximation. 
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I. INTRODUCTION 



In recent years, much attention has been attracted by 
the so-called short-time critical dynamics iSiSi^iSiSiLSiS A 
thermodynamic system in a critical point is characterized 
by long-range correlations in equilibrium, so when such 
a system is quenched from a high temperature Tq ^ 
to the critical point Tc, the growth of correlations gov- 
erns the relaxation process. As shown in Ref. Q, the 
relaxation process displays a universal scaling behavior 
already at times just after a microscopic time scale tmic 
needed for the system to forget its microscopic details 
and lasts until the system has lost all memory about its 
macroscopic initial condition. The initial nonequilibrium 
stage of the critical relaxation is characterized by new 
critical exponents 6 and 9' , which describe the behavior 
of the response function G(t, t') oc {t/t'Y for t' ^ and 
the initial increase of the order parameter m{t) oc , 
respectively. The critical exponents 9 and 9' depend on 
the dynamic universality class ,JA and have been calcu- 
lated for a number of dynamic models such as the model 
with a nonconserved order parameter^ (model A), the 
model with an order parameter coupled to a conserved 
density^ (model C), and the models with reversible mode 
coupling'^ (models E, F, G, and J). The universal scaling 
behavior of the initial stage of the critical relaxation has 
been verified by extensive numerical simulations4i 

The influence of various kinds of disorder on phase 
transitions is one of the central problems in condensed- 
matter physicsiii*i2iiiiiii^i^ In the present paper, we 
focus our attention on the nonequilibrium critical prop- 
erties of systems in the presence of quenched temper- 
aturelike disorder (for example, magnetic systems with 
quenched nonmagnetic impurities). In the simplest case, 
the disorder may be viewed as randomly distributed 
pointlike defects>i& According to the Harris criterionjii 
the quenched uncorrelated pointlike defects change the 
critical behavior only if the heat capacity exponent of 
the pure system is positive ap > 0. As a result, the 
pointlikc disorder is relevant only for Ising-like three- 
dimensional (d = 3) systems. The short-time critical 
behavior of systems with pointlike uncorrelated defects 
was studied in Refs. and |^ Real systems, however, 



often contain defects in the form of linear dislocation, 
planar grain boundaries, three-dimensional cavities, or 
other extended defects. One possibility to treat the sys- 
tems with extended defects is the model suggested by 
Weinrib and Halperin (WH) ^ in which defects are long- 
range correlated and characterized by a correlation func- 
tion that has a power-law decay g(x) cx x~°' with distance 
X. This type of disorder has a direct interpretation for 
integer values a: the case a — d corresponds to uncor- 
related pointlike defects, while a = d— l{a = d— 2) 
describes infinite lines (planes) of defects of random ori- 
entation. The statics of the WH model and its dynamics 
not far from equilibrium were examined by means of the 
renormalization-group (RG) methods to two-loop order 
using both the double e — A — d, (5 = 4— a expansior*iS 
and the direct calculations in d = 3*^ The nonequilib- 
rium critical relaxation of the WH model was considered 
only in the one-loop approximation,^ although numerous 
investigations of pure and disordered systems performed 
with the use of the field-theoretic approach show that 
the predictions made in the one-loop approximation, es- 
pecially on the basis of the e expansion, can differ from 
the real critical behavior liSi^S The possible alternative to 
the isotropic scenario realized in the WH model is the 
anisotropic scenario realized in the model proposed by 
Dorogovtsev^ in which defects are strongly correlated 
in Sd dimensions and randomly distributed over the re- 
maining d = d — Ed dimensions. The case Ed = is 
associated with uncorrelated pointlike defects, while the 
extended linear (planar) defects are related to the cases 
Ed — 1(2). The case of the noninteger value of £d may 
be related to a system containing fractal-like defects so 
that Ed is interpreted as an effective fractal dimension 
of a complex random defect systemi2^ The statics and 
equilibrium dynamics of Dorogovtsev's model were stud- 
ied in Refs. 23-24 25-26 (for a review of the models with 
extended defects, see Refs.iHilli). 

In the present paper, we investigate the short-time 
critical dynamics of Dorogovtsev's model with a noncon- 
served order parameter using a field-theoretic approach 
in the two-loop approximation. The paper is organized as 
follows. Section introduces the model describing the 
critical dynamics of disordered systems with extended 
defects. In Sec. IIIII the corresponding effective field the- 
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ory is renormalized up to two-loop order. In Sec. II VI we 
derive the asymptotic behavior of response and correla- 
tion functions and calculate the critical exponents using 
the Pade-Borel resummation method. The final section 
contains our conclusions. 



II. THE MODEL 

In equilibrium at temperature T the Hamiltonian de- 
scribing the disordered systems with extended defects is 
given by2i 



Hi 



1 " 



/3=1 



IV 



p,2, _^ 9^ 



13=1 



,(1) 



where is the n-component order parameter, r is 

the reduced temperature, and g is a positive constant. 
a is the bare anisotropy constant and the factor a^"* in 
the quartic term is included for technical convenience^ 
V{x.±) is the potential of defects, which can be viewed 
as e^-dimensional objects, each extending throughout the 
whole system along the coordinate X|| , while in the trans- 
verse direction x^ they are randomly distributed with 
the concentration taken to be well below the percolation 
limit. The potential V{x±) is assumed to be Gaussian- 
distributed with the zero mean and the second cumulant 



(F(xi)F(xl)) = A<5'^-^''(xi - x'^ 



(2) 



The dynamics of the non-conserved order parameter 
4>^{x) in the system defined by Hamiltonian ^ can be 
expressed in the form of the Langevin equatioufiS 



dt 



where A is the Onsager kinetic coefficient, and the func- 
tion rj^ (x, t) is a Gaussian random- noise source with cor- 
relations 



7?'^(x,t)7?^(x',f) = 2A^/3-y(5(x-x')(5(t-i')- (4) 

Here, we summarize the main features of model 
revealed in Refs. and I2I2I25II2I The Harris crite- 
rion is modified in the presence of extended defectsi^^ 
The disorder affects the critical behavior only if the cor- 
responding crossover exponent Lp is positive, 



if = ap + Vp£d > 0, 



(5) 



where Vp is the correlation length exponent of the pure 
system. Consequently, the disorder with extended de- 
fects is relevant for d = 3 over a wider range of n than the 



pointlike disorder.-— The critical properties of the model 
can be examined within the RG framework by using the 
double exp ansion in both e and Sd, which was suggested 
in Ref . Due to the spatial anisotropy caused by disor- 
der, two correlation lengths and ^|| naturally arise, one 
of which is perpendicular to the extended defects direc- 
tion whereas another is parallel to this direction. In the 
critical point, their divergences are characterized by cor- 
responding critical exponents and : S^±_ cx |t|^''^, 
^11 cx |t|~''ii . The correlation of the order-parameter fluc- 
tuations in two different points depends on the orienta- 
tion of their distance vector, so that the behavior of the 
correlation function is characterized by a pair of Fisher 
exponents r^^ and Tyy. As was shown in Ref. |24L the equi- 
librium dynamics is also modified and there are two dy- 
namic exponents z_l and z\\ . On the other hand, as the 
interaction of all order-parameter components with dis- 
order is the same, the susceptibility (as well as the order 
parameter and heat capacity) is characterized by the sin- 
gle exponent 7 and can be written a32i2iSSiS& 

x{q^,q\\,t,T)^T-<g, (r-"-- , r'^^^ g,i , tr^^ ) , (6) 

or in the critical point (r 0) as 

X(gi,9||,t,0) = \qx_\'^-^Q2{q\\\qA-''^l''\t\q^Y-) 

- kiir'"-^g3(g±kiir^"/'^\tkiir"), (7) 

where q^, q\^ are the components of the momenta along 
the x^ and X|| directions, respectively, and Qi are the 
scaling functions. The critical exponents defined above 
were calculated to second order in e and Ed in Refs. 12^ 
andlH 

It has been recently argued in Ref. that the planar 
{sd = 2) defects can destroy the sharp continuous tran- 
sition in the Ising-like systems due to the existence of 
rare infinite spatial regions which are devoid of defects 
and therefore may be locally in the ordered phase. It 
had been proposed early that the similar effects can give 
rise to the instability of the usual critical behavior of 
disordered systems with respect to the replica symmetry 
breaking^ (RSB) and even lead to the appearance of an 
intermediate spin-glass phase. '^■^ However, for pointlike 
defects the rare regions have finite size and thus cannot 
develop the true static order, so that the order parameter 
on such rare regions still fluctuates. The more accurate 
RG calculations performed in Ref. have shown the 
stability of the critical behavior of the weakly disordered 
systems with respect to RSB potentials. In the case of 
extended defects, these regions are infinite in £d dimen- 
sions and, according to the Mermin- Wagner theorem^M 
can develop the true static long-range order for n — 1 and 
ed > 2 [in two-dimensional systems with continuous sym- 
metry (n > 1) there is no true static long-range order]. 
As suggested in Ref. [s^, the last can result in the smear- 
ing of the phase transition. Although this picture, based 
on the extremal statistics, is supported by the accom- 
panying numerical simulations, "^"^ there is a contradiction 
with the early numerical studies. '^^ Considering in what 



follows the particular case n = 1 and = 2, we will not 
take into account these effects. 

We now consider the relaxation of model Q-Ql from 
a nonequilibrium state with a small correlation length. 
This initial state can be macroscopically prepared by 
quenching the system from a high temperature Tq 3> Tc 
to the critical temperature Tc- We have to specify the dis- 
tribution of the initial condition 0o(x) — (j)^{x,t — 0). 
The reasonable assumption is that the distribution is 
Gaussian, 



cxp (^-^ J d'xi^^.i^) 



<(x))^ 



(8) 



where mQ(x) is the initial order parameter and Tq^ is 
the width of the initial distribution. Equation ^ guar- 
antees that the initial correlations are short-range. Terms 
of higher order in the exponent of Eq. (jS)) turn out to be 
irrelevant in RG sense. By introducing a Martin-Siggia- 
Rose response field^^ (f)^{x,t), we can average over the 
random noise ri^(x,t) and obtain an equivalent formula- 
tion of dynamics in terms of a generating functional for 
all response and correlation functions, 



W\Kh] = 



dt [ d'^xih^cj)'^ + hf^^^) 



(9) 



where the action Cv[4>, 4>\ is given by 



dt / d'^x4>'^ 



dt 



(10) 

We can perform the average over disorder directly in 
Eq. lO without using the replica trick^^ and obtain the 
T^-independent and translational-invariant (however, in 
space x only) effective action 



C[(j),$\ = dt d'^x4>'^ 









FIG. 1: Diagrams contributing to To,! (g, at two-loop or- 
der. Diagrams Hi — Ha contribute to Tq (g, ; diagrams 
Ri,R2, and R3 contribute to K{q; t' , t"). The solid lines corre- 
spond to the initial correlator c'^\t,t') in diagrams H\ — He 

and to the equilibrium correlator Cq'^^ {t, t') in R\,R2, and Rz; 
the lines with an arrow correspond to the propagator Gq{t, t'); 
the four-leg vertex x corresponds to the g interaction; the 
wavy line corresponds to the disorder interaction A; the wall 
corresponds to the "time surface" t — 0. 



then given by 

Gq{t, t') = Q{t - t') exp(-A(T + a^ql + q\){t - t')){12) 
C,{t,t') = C^^\t,t') + Cj;\t,t'), (13) 
where 

exp{~\{T + a^q?.+ql)\t-t'\) 



T + a^q}, + q\ 



Cfit,t') = (r^-i- 



r -I- a^q?. + q\ 



X exp(-A(T + a2q2 +g2^)(t + t')) (15) 

are the equilibrium and the initial nonequilibrium parts 
of the Gaussian correlator, respectively. 



III. RENORMALIZATION 



- vi) / + \^^'^<t>'^^'' - X4>^ - Ix'A Jdt' 



X / dVj''(xi -x'^)/(x,t)0'^(x',t')</)^(x',i') 



(11) 



For g = A = 0, the generating functional ^ becomes 
Gaussian and can be easily evaluated in momentum space 
by solving the corresponding variational equations. We 
have to take into account the initial condition JHJ by 
imposing the boundary conditions 0''(x, 00) = and 
0q(x) = TOq (x)-f rQ'"'^0''(x, 0). The free response function 
(propagator) Gq{t, t') and the free correlator Cq{t, t') are 



We now exploit the methods of renormalized field the- 
ory in conjunction with a generalized expansion in e to in- 
vestigate the scaling behavior of nonequilibrium response 
and correlation functions. We define the full one-particle 
reducible Green functions as 



G 



(16) 



The Green functions ^ may be calculated by a per- 
turbation expansion in the coupling constants g and A. 
The Feynman diagrams that contribute to the Green 
functions computed using the effective action (|ll|l in- 
volve momentum integrations of dimensions d = 4 — e 
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TABLE I: Values of two-loop diagrams contributing to 
ro.i(g = 0, The common weight factor for all diagrams 

is (n + 2)/12{Kag){K^A)\r{l - (e + e)/2)(2At)-i+(^+^"'/2. 



Hi 


— 1/e — 3/4 — £/4£ 


H2 


l/£ 


H3 


l/£ - ln2/2 


Hi 


2/£ + 1 - 21n2 


H5 


2/£ + 1 + 51n2 - 31n3 - \/31n(2 + V3) 


He 


-l/£ + ln2/2 


Ri 


2/£ 


R2 


2 In 2 


R3 


3 In 3 - 4 In 2 



and d = 4— e. We use the dimensional regularization- to 
calculate integrals and employ the minimal subtraction 
scheme^ito absorb the remaining poles in (c\£~^C2s) into 
multiplicative Z factors and to introduce the renormal- 
ized quantities according to 

Tb = (Z,/Z,)r, Afc = (Z,/Z,)1/2a, afc = Z„a,(17) 

where the subscript 'b' denotes bare quantities and fi is 
an external momentum scale. Naive dimensional analysis 
gives To oc n"^ and, thus, Tq^ is an irrelevant parameter 
in the RG sense. It is convenient to consider the Dirich- 
let boundary conditions Tq^=0 and mo(x) = 0. The 
general case is recovered by treating the parameters Tq^ 
and r7io(x) as additional perturbations. The Z factors, 
except for Zq, were calculated in Ref . [25I to two- loop or- 
der. The new factor Zq serves to cancel the divergences 
arising from the initial part Cq'"* (t, t') for t + t' —>■ 0. Note 
that although there exist two different Z factors for fields 
4> and 0, we need only Zq to renormalize fields (j^o and 
(j>o- This is due to the Ward identities 

0^(x)=ro-i0(?(x), (18) 
9t/(x,t)|t=o = 0[^(x) = 2A0(^(x), (19) 

which hold when inserted in the connected Green 
functions.^ In order to determine Zq, we require that 

(ZoZsZ^-)"^/2Gj,i(x,t) = finite fore, e ^ 0. (20) 

It is more convenient to work with the Fourier trans- 
form Gl i{q,t) — J d'^x exp{—iqx.)Gl i{x,t) , which can 
be written in the following form: 

t 

GlAq,t)^ j dt'G^.,{q-t,t')T'i^^,{q,t'\^^^, (21) 


where Gi^i{q]t,t') denotes the full unrenormalized two- 
point function coming only from the equilibrium part of 



the correlator ifTH) . ro'\(9, t')[<^o] ^® ^ reducible one-point 
vertex function with a single insertion of the response 
field 00 5 which contains in its last irreducible part at 
least one initial correlator H15|l . Figure 1 shows the two- 
loop diagrams Hi ~ Hq which contribute to To ^')[<ij 
together with those computed for pure systems in Ref.|2 
The singular parts of computed diagrams are given in 
Tabled Note that although Gi^i{q;t,t') is calculated 
with equilibrium propagators and correlators, it is differ- 
ent from the translational-invariant equilibrium response 
function G^i^^ {q;t — t') because of the restriction of inter- 
nal time integration to positive times t > 0. Indeed, the 
equilibrium Green functions must be calculated by means 
of the effective action (lllll with the lower limit of time 
integration t — replaced by i = —00. We can relate 
Gi,i(g;i,t')to G^^f{q,t-t') using the method suggested 
in Ref.0- Following Ref.|3, we choose the equilibrium dis- 
tribution Poq[</'o] exp(— Hy[0o]) instead of distribution 
(|H1) to average over the initial field 00 in Eq. (jSJ. This al- 
lows us to obtain a perturbation series for G^if' (g, t — t') 
where internal time integrations range from zero to infin- 
ity. However, the non-Gaussian probability distribution 
Pcq [0o] generates in the effective action (|11() an additional 
vertex, 

- ^AA j dt j d'^x j dVj'^~(xx -x'_l) 


X 4>''Mcp('{^,t)4>l{^)4PQ{^), (22) 

which is located at "time surface" t — Q. We can re- 
late G'^^f (g, t-t')io Gi,i ((j; t') by using the following 
equation: 

t 

Gt.f{q,t~t')^ J dt''GiA<i;tXKfiqX\^^,)^, 

(23) 

where T^f'^q, ^")[0(t')] is a reducible vertex function with 
a single insertion of (j){t'), which is calculated with equi- 
librium correlators and contains in its last irreducible 
part at least one of the new vertex H22|l . In contrast to 
the pure system,^ the disordered system has al- 

ready the nontrivial contributions to ro°i''(9' ^")[0(t')] 
two-loop level. The corresponding two-loop diagrams Ri 
- i?3 are given in Fig. 1. Considering Eq. H23|l as an inte- 
gral transformation with kernel ro°i''(5, ^")[0(t')]i can 
find the kernel of the inverse transformation K{q] t, t') 
order by order in perturbation theory using the following 
integral equation: 

t 

J dt'rifiq,t\^^^,,^^Kiq-X,t')^S{t-t'). (24) 
t' 
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Combining Eqs. (|?T|l . and (|^ . we arrive at 
t 

Gj,i(g,i) = 1 dt'Gl';?^(<z,t-t')ro,i('Z,0[^o]' (25) 



where 

t' 

ro.!!?,*')!^,] = / dt"if(g;t',t")4:\(g,t")[^„]. (26) 



The two-loop result for the singular part of vertex func- 
tion (|26|) is given by 



ToAQ = 0,tr^^^ - <5(t) + ^^Agr (l - 

- ^^Ag^r^ f 1 - I 
36e V 2 

(n + 2) 



{2Xt) 



-l+s/2 



r(l+e) "2,' ^2^^) 



+ 



^2 Ag^r(l-£) (^-- + ln2--) (2At) 



-1+e 



(n + 2) 



48 ^ V 2 



24 



A-2- 



(2Ai) 



-l + (£+e)/2 



(27) 



r 



where A = 7 -h 41n2 - 4\/31n(2 V3), r( x) is the Eulcr 
gamma function, and we have absorbed factors of = 
27r^/V((27r)''r((i/2)) and into the redefinition of the 
coupling constants g and A, respectively. Taking into 
account Eq. H25() , the renormalization condition H20|l can 
be rewritten in the formii^ 



-1/2 



dte"'"*ro,i(o,t)^|;s 



finite for e, e ^ 0, 



(28) 

where Fq,! (0, is given by Eq. l(77|) expressed in terms 

of the renormalized quantities using Eqs. p7ll . They read 
in explicit form 



9 ■ 
A 



5 1 



A/, = A 1 



in + 8)g 6A , ~ ~2n 

6e £ 



h 0(e^,ee, e^) 



As a result, at two-loop order we obtain 



Zn = 1 



6e 



12e2 



ln2--)s 



2 (e2 - 24e + Aee) 



24 



ee(e -I- e) 



5A. 



(29) 



For n = I and e = e, Hamiltonian J^l corresponds to 
the Ising model with quenched pointlike disorder, so that 
Eq. H29|) reduces to Zq obtained in Ref. |^ 



IV. SCALING AND CRITICAL EXPONENTS 

We are now in a position to discuss the scaling proper- 
ties of Green functions p6|l . Using the fact that the bare 
Green function G^x does not depend on the external 

N,Nb 

momentum scale fj, introduced in Eqs. (|17|l . we can derive 
the RG equation in the usual way.'^^ It reads 



. N N L ^ 
+ y7 + y7 + -^(7 + 7o) 



G 



N,N 



0, (30) 



where we have introduced the operator R = fidf^ + Pg^g + 
f^AdA+lTTdr+jxXdx+jaCtda, and the P and 7 functions 
are defined as derivatives at constant bare parameters, 

Pg = ^J■^^, Inglg , (3a = f^df. In A|g , 7^ = Inal^ , 
7A = A*9^1nA|o = (7 - 7)/2, 7^=^9^1nT|o, (31) 
7 = ^idf. In Z^Iq , 7 = fidf. In Z^l^ , 70 = ndf. In ZqIq . 

Functions H31(l , except for 70 , were calculated in Ref. 
to two-loop order. The new function 70 reads 



70 



6 



1 



In 2 - - 
2 



4 V e 



(32) 

The nature of the critical behavior is determined by the 
existence of a stable fixed point (FP) satisfying the si- 
multaneous equations 



Pg{g\A*) = 0, /?A(ff*,A*)=0. 



(33) 
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The full set of F Ps exhibi ted by Eqs. was discussed 
in Refs. 113 and I23l24l25i We are interested in the FP 
which controls the critical behavior affected by extended 
defects. The corresponding FP is given for n > 1 bjiS^ 



9 



3(3e- 2e) 



2(n-l) 256 
+4(73^2 + 52n + 352)ee - (llln^ 
-32(n + 2)(n 

{n + 8)e-2{n + 2)s 



-4(3l7i2-68n+196)£^ 

696n + 624)e2 



l)e3e-i)(n-l)-3, 



A* = 



(n-1) 



(34) 



In the case n = 1, there is an accidental degeneracy in 
Eqs. ijSSI), which leads to FPs of order e^/^ and e^/^. The 
FP corresponding to the systems with extended defects 
reads 



9 



4A* = 8[(9e-6e)/106; 



1/2 



(35) 



The solution of Eq. H30|l and the simple dimensional 
analysis yield the scaling behavior of the Green functions 
at FPs, 

g|^^({x|| , x^, 0, r, A, a, g*, A*; ^) = /e.(i-C)(^+^+L-2)/2 

^ l(d-2+rtj_)N/2+{d+2+fi)N/2+{d+2+7io+fi)L/2 

X G|^^({x||/Sx^/,tP-},r/-i/--,A,a,5*,A*;/.), (36) 
where the critical exponents are given by 

C=l-7a(ff*,A*), 1/^.^ = 2-7,(3*, A*), 

77^ =7(g*,A*), = 2 + 7^,(3*, A*), (37) 

V^j{g*,A*), % = 7o(5*,A*). 

The first factor on the right hand side of Eq. H36|) arises 
from the dependence of Green functions on a not only 
through X|| , but also through an overall factor, which 
can be easily found out by inspection of the Feyn- 
man diagrams (see the Appendix of Ref. [25j) . Using 
Eq. (|36|) . we reveal the known relations for the critical 



21.23.24 



exponents^ 

= (2-77^)/(2-?7||) =zj_/z\\ 



(38) 



Substituting FP into Eq. (|^ . we obtain the new 
critical exponent 770 for n > 1 as a double expansion in e 
and e. 



r]a = (n + 2)(n- l)"M0-5e-0.75e + (n- 1)- 

X ((0.201599^2 + 0.601119n + 1.992204)e2 

- (0.553223^2 - 0.400168n + 3.573508)ee 

+ (0.284424^2 - 0.684984n + 1.642748)e2 

+ (0.03125^2 - 0.0625 + 0.03125)e^£~^)] . 



In the case n = 1, we have 

r/o = -4[(9£- 6e)/106]i/2^ 



(39) 



(40) 



Note that Eq. (|39|l involves the ratio e/e, which is of order 
unity within our approximation. The authors of Ref. [3 
pointed out that this ratio, which appears in g* and 7, , is 
absent from the final critical exponents z.j_, ri±, z±, and 
C at second order. This encouraged them to believe that 
similar cancellations will occur at all orders. It is easy to 
see that the reason for the ratio cancellation is that the 
last three exponents do not contain the first-order term 
proportional to g* , which gives rise to the ratio e/e at sec- 
ond order. In the case of vx, the accidental cancellation 
of terms e/e occurs when one substitutes FP into 7,. 
However, the cancellation can be destroyed, for example, 
by introducing cubic anisotropy^Sfi so that all exponents 
involve the ratio e/e. For physical applications, we are 
interested in positive values e and e = e + ed, for which 
< e/e < 1, and therefore the critical exponent ifll^ is 
well defined. 

In order to describe the scaling behavior of the 
nonequilibrium critical relaxation, we employ a short- 
distance expansion of the fields </> and </> in terms of the 
initial fields^ for t ^ 0, 



(41) 



Introducing the rcnormalized amplitude functions ac- 
cording to ab{t) ~ ^^^cr(i) and dh(t) = Z'^^^'^ait), 
we can derive the corresponding RG equations [i? — 
?7o/2]cr(t) = and {R - ^Q/2\a{t) = 0. Taking into ac- 
count [a] ^ {Xfi^)~^ and [a] ~ 1, we find the solutions of 
the RG equations at FPs 

a{t, T, A, a; n) = l-'^-""/^cr{tF^ , rr^/"^ , A, a; fi), 

a{t,T,X,a;n) = l-""^^a{tF^ ,tI-^^''^ , X,a; fi). (42) 

We are interested in the short-time critical behavior of 
the response and correlation functions. Inserting ex- 
pansions (|41|) into appropriate Green functions, we ob- 
tain G{x,t,t') = CT(t')Gj i(x,t) -f ... and C{x,t,t') = 
2Xa{t')Gl i{x.,t) + where we have used Eq. |[T!?|l . 
Combining Eqs. (|36|l and 1)42(1 . we obtain the scaling 
behavior of the response and correlation functions for 
^ and ^I'^t > 0, 



C(x,t,i') = 4"''""^ 



1 1 p ( ^ ^ J— 

t') xi'CW 



,(43) 



where the new dynamic critical exponent is defined by 

e = -Vo/{2z^)- (44) 
To second order in e and e, it reads 
e = (n + 2)(n- l)-i [0.1875e-0.125£- (n- l)-2 
X ((0.062119n2 -|- 0.232311n + 0.404301)e2 
- (0.169556^2 - 0.021917n + 0.784002)ec 
-I- (0.007813n2 -|- 0.007813n - 0.015625)e^e"^ 
-I- (0.086731^2 - 0.155621n + 0.379437)e2)] . (-45^ 
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TABLE II: Values of critical exponents rjo, 6, and 9' for d — 3. 





^ a 


^70 


y 


a 


\ 


0" 


—0.374 


0.092 


0.087 




1 


—1.346 


0.252 


0.140 




2 


— 1.780 


0.308 


0.331 


2 


0' 


-0.647 


0.162 (0.160) 


0.200 




1 


-0.614 


0.142 (0.138) 


0.059 




2 


-0.688 


0.157 (0.152) 


0.054 


3 


0' 


-0.790 


0.198 (0.196) 


0.229 




1 


-0.834 


0.191 (0.186) 


0.111 




2 


-1.008 


0.226 (0.214) 


0.111 



"Taken from Ref. 

''Computed using the results of Ref. (Ij for pure systems. 



Although the system is characterized by the single expo- 
nent 6, the scaling laws (|43|l reflect the strong anisotropy 
of the system under consideration. 

We now discuss the scaling behavior of the order pa- 
rameter to'^(x, t) as a function of its initial value mg (x). 
For the sake of clarity, we will suppress the superscript 
f3 in what follows. Considering mo(x) as an additional 
time- independent source coupled to the fleld (/'o(x), we 
add the term J c?'^a:;mo(x)(^o(x) to the effective action 
(|ll|l . The order parameter is given then by 

to(x, t) = (0(x,i)) = 5W[h, h]/6hi^, t)\^=h=o- (46) 

Expanding Eq. H46|l in ma we obtain for the case of a 
homogeneous initial order parameter 

m(mo,t) = X]Ty / d'^xi... f d'^XLG^^,{{x},t)m^ . 

L=l ■ 

(47) 

Equation (|47|l holds also for the renormalized quantities, 
so that substituting Eq. (|36|l into Eq. (|47l) and summing 
over L, we obtain 

m{mo,t) = mot"' F{mot'''+^^/''^'^,Tt^/''^'^), (48) 

where 9' = -[r]j_ + f] + rjo + 26^(1 - C)]/{2z±), (3 is the 
critical exponent for the order parameter and we have 
used the scaling relation 2f3 = i^±{ri± — 2 + rf) + edi'\\-'^^ 
For n > 1, we have 

e' = 0.125£-|- (n- i)-2 

X [(0.058228n2 - 0.062655n + 0.204281)£e 

-t- (0.024597n2 + 0.055126n + 0.011862)£2 

- (0.007813^2 -I- 0.03125n + 0mi2b)e^r^ 

- (0.034287^2 -t- 0.088644n + 0.145612)e2] . (49) 

The critical exponents 6 and 9' are related by 

e' ^0 + i2-T^^-z^-Sd{l-O)/z±- (50) 



The scaling function F{yi,y2) has a universal behavior 
at the critical point r = 0: -F(0, 0) is finite, while for 
Ui oo, F{yi,0) cx 1/yi. Thus, Eq. exhibits the 
crossover between the initial and the equilibrium stages 
of the critical relaxation at time tc oc TOq ^^'•^ +l3/iy±z±) ^ 
when the initial increasing m{t) cx t" turns to the long- 
time decreasing m{t) oc t~^/'^^^^. 



It is well known that the convergence of double series 
in £ and e is poor.^^'^e _ \^^fJQ know only the lowest 

order term 1)4011 and cannot apply any resummation tech- 
nique. Therefore, we compute exponent (|40|l as it is, and 
then we calculate and 9' by using Eqs. (|44)l and 1)50)1 . 
In order to estimate the values of the critical exponents 
770, 0, and 6' for n > 1, we employ the Pade-Borel resum- 
mation method extended to the two-parameter case (see 
Ref. 15 and references therein). Numerical values of the 
critical exponents 770, 0, and 9' computed for the three- 
dimensional Ising (n — 1), XY [n = 2), and Heisenberg 
[n = 3) systems with pointlike (e^ — 0), linear (e^ = 1), 
and planar (e^ = 2) defects are presented in Table ^ 
To calculate 9 and 6*', we can use either expansions 1)45)) 
and 1)49)1 . or scaling relations 1)44)1 and ()5())l . This can be 
utilized to check the accuracy of the resummation proce- 
dure. The values of 9 given in parenthesis (see Table lll|l 
are computed with the use of Eq. ()44)l . The small discrep- 
ancies indicate the fairly good accuracy of the computed 
values of 9. Unfortunately, using the scaling relation 1)50)1 
leads to unreasonable values of 9' . The reason is the large 
error in the determination of rj±. Indeed, the direct cal- 
culation of rj±_ gives 77^ w —0.44 for n = 2 and = Ir^ 
while the scaling relation 77^ = 2 — 7/1^^, where 7 is the 
critical exponent of susceptibility, suggests 77^ ~ 0.008ii2& 
The resummation of expansion ()49)l results in more plau- 
sible values of 0' , which are shown in Table Hll 



According to the Harris criterion^ the pointlike disor- 
der {sd = 0) is relevant for the critical behavior of three- 
dimensional systems only if n = 1. The extended defects 
change the critical behavior for all n li Sd > £^^"^^(71), 
where the marginal values £^^"^^(1) — —0.173, £^^"^^(2) = 
0.016, and £^^"^^(3) = 0.172 were calculated for d = 3 in 
Ref. |2y using the extended Harris criterion Q • As one 
can see from Table IhI the 3D Ising system with extended 
defects is characterized by larger values of the initial slip 
exponents 9 and 9' than the 3D Ising system with point- 
like disorder. For ti > 1, the pointlike defects do not af- 
fect the critical behavior and the critical exponents take 
the same values as that of the pure system. Table IhI 
shows that for ti = 2, 3, the presence of extended defects 
leads to the small changes in 0, while the corresponding 
decrease of 9' is more significant and can be observed in 
numerical simulations. Remarkably, the decrease of 9' 
caused by extended defects almost does not depend on 
Ed for n > 1 in the considered approximation. 
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V. CONCLUSIONS 

In the present paper, we have studied the short-time 
critical dynamics of the system with extended defects, 
which are strongly correlated in Ed dimensions and ran- 
domly distributed over the remaining d = d — Ed dimen- 
sions. We have considered the nonequilibrium relaxation 
of the system starting from a macroscopically prepared 
initial state with short-range correlations and found the 
scaling behavior of the nonequilibrium correlation and 
response functions, and the initial increase of the order 
parameter. Using the Pade-Borel resummation method, 
we have calculated the initial slip exponents 9 and 0', 
which describe the initial growth of correlations, to two- 
loop order. 

Unfortunately, to our knowledge there is only one nu- 
merical investigation of the short-time critical dynamics 
of systems with extended defects. In Ref. |3^ the short- 
time critical dynamics of the 2D Ising model with bond 
disorder perfectly correlated in one dimension and un- 
correlated in the other was studied using Monte Carlo 
simulation. This model corresponds to the case e = 2 
and Ed = 1, for which Eq. ^ gives 9' = 0.154. The 
value 9' — 0.176(2) reported in Ref. |^ is larger than 
the prediction of the RG analysis. The difference can be 
attributed to very strong disorder used in Ref. [s^ (the 
concentration of bonds p = 0.5), while the RG descrip- 
tion is valid only for the regime of weak disorder (p closes 
to 1), in which the asymptotic behavior does not depend 
on the concentration p.^^ 

The relevant question is how the critical properties of 
the system are modified if the defects, oriented along the 
direction xy , are very long but finite, instead of being 
extended throughout the system. According to the usual 
scaling theory near the critical point Tc, the only relevant 
scale that remains in the system is the correlation length, 
which diverges in the critical point as ^ oc |t|^'^''. If the 
disorder is weak, its effect on the critical behavior in the 
vicinity of the critical point is negligible as long as the 
correlation length ^ is smaller than the average distance 
between the defects Zdof, i-e., the critical behavior is con- 



trolled by the FP of the pure system. In close vicinity of 
the critical point (t 0), the correlation length grows 
and becomes large than ^dcf for t ^ JSr'^/v ^ where Lp is 
the crossover exponent ©. However, it is possible that 
the correlation length remains smaller than the length of 
defects. For these values of r, the critical behavior is con- 
trolled by the FP of a system with extended defects, the 
two correlation lengths naturally arise due to the strong 
anisotropy, and one can use the results obtained here for 
the case of the strong correlated disorder. Closer to the 
critical point, the correlation length becomes larger 
than the length of defects and the asymptotic behavior 
exhibits crossover to the critical behavior of the system 
with pointlike disorder. The corresponding crossover ex- 
ponent is I'll . 

From an experimental point of view, the considered 
model may be relevant for the anomalous scaling be- 
havior of the so-called narrow component in magnetic 
systems such as Ho and Tb, which exhibit two different 
length scales for critical fluctuations usually ascribed to 
long-range correlated quenched disorder.'*'^ Other possi- 
ble applications are high-Tc superconducting materials, 
where extended disorder is caused by irregularly alter- 
nating layers, such as those containing Y and Ba in Y- 
Ba-Cu-0 structures, or CuO chainsiSI The obtained re- 
sults can be used to investigate the long-time dynamics 
of systems with extended defects in the aging regime^i 
They can be helpful also for numerical simulations of sys- 
tems with extended defects, especially by using new effec- 
tive algorithms, based on the fact that due to the small 
correlation length during the initial stage of the critical 
relaxation, the determination of critical exponents re- 
quires less effort than numerical studies of equilibrium 
dynamics4^ 
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